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We review recent progress in high-resolution imaging of scattered light from disks around 
young stellar objects. Many new disks have been discovered or imaged in scattered light, and 
improved instrumentation and observing techniques have led to better disk images at optical, 
near-infrared, and thermal-infrared wavelengths. Multi-wavelength datasets are particularly 
valuable, as dust particle properties have wavelength dependencies. Modeling the changes in 
scattered-light images with wavelength gives direct information on the dust properties. This has 
now been done for several different disks. The results indicate that modest grain growth has 
taken place in some of these systems. Scattered-light images also provide useful constraints on 
the disk structure, especially when combined with long-wavelength SEDs. There are tentative 
suggestions in some disks that the dust may have begun to settle. The next few years should see 
this work extended to many more disks; this will clarify our understanding of the evolution of 
protoplanetary dust and disks. 



1. INTRODUCTION 

This chapter reviews the progress since the Protostars 
and Planets IV meeting (McCaughrean et al., 2000) in ob- 
servations and modeling of scattered-light images of disks 
around young stellar objects (YSOs). For many years these 
disks were studied in the infrared and at millimeter wave- 
lengths, but without the spatial resolution necessary to re- 
veal their detailed structure. With the advent of the Hubble 
Space Telescope (HST) and ground-based adaptive optics 
(AO) systems, these disks were seen in scattered optical and 
near-infrared starlight. The first images of disks in scattered 
light showed them as never before. These images currently 
provide the highest spatial resolution images of disks and 
provide unique information on disk structure and dust prop- 
erties. 

Most disks around young low- and intermediate-mass 
stars fall into one of two categories on the basis of their 
gross observational properties. YSO disks are optically 
thick at visible and near-infrared wavelengths, are rich in 
molecular gas, and are found around Class 1 and Class 11 
systems. Debris disks are optically thin at optical and near- 
infrared wavelengths, have only trace quantities of gas, and 
are found around Class 111 and older systems. 

Our current understanding of the evolution of dust in 



disks around young low- and intermediate-mass stars is that 
these observational categories correspond to very distinct 
phases. In YSO disks our understanding is that gas and 
dust from molecular cores is processed through the disk and 
provides the raw material both for accretion and outflows. 
The dust is processed first in the dense molecular core and 
then in the disk itself. The dust grows and suffers chemi- 
cal processing. Dust growth is hypothesized to result in the 
production of planetesimals, which form rocky planets and 
the cores of gas giant planets. In debris disks, our under- 
standing is that dust is present largely as a result of colli- 
sions between planetesimals. Thus, YSO disks are thought 
to be characterized by dust growth whereas debris disks are 
thought to characterized by planetesimal destruction. 

In this chapter we focus on observations and modeling 
of scattered-light images of YSO disks around low- and 
intermediate-mass stars. That is, we focus on the phase in 
which dust is expected to grow. Of course, scattered-light 
observations and models are not the only means to explore 
these disks. The chapter by Dutrey et al. discusses mil- 
limeter and sub-millimeter observations of the gas and dust 
components and the chapters by Bouvier et al. , Najita et al. , 
Millan-Gabet et al. cover observations of the inner disk re- 
gion, notably interferometry, and their interpretations. The 
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Fig. 1. — Scattered light images of four face-on or intermediate inclination YSO disks newly resolved since 1998. Left-to-right: TW 
Hya face-on disk (Krist et al, 2000), HD 100546 (Grady et al, 2001), HD 141569A (Clampin et al, 2003), and HD 163296 (Grady et 
al, 2000). 



chapter by Monin et al. discusses aspects of disks that are 
peculiar to binaries. Further from the subject of our chapter, 
the chapter by Cesaroni et al. discusses disks around young 
high-mass stars and the chapter by Meyer et al. reviews de- 
bris disks around solar-type stars. 

In the following sections we summarize the data avail- 
able on scattered-light disks around YSOs, imaging tech- 
niques, and what we can learn about disk structure and 
dust properties from modeling multi-wavelength images 
and SEDs. We also speculate about future advances. 

2. OBSERVATIONAL PROGRESS 

At the time of the Protostars and Planets IV meeting, 
held in 1998, only a dozen YSO disks had been resolved in 
scattered light and another half dozen in silhouette against 
the Orion nebula. Since then, the number of resolved disks 
in these two categories has doubled and quadrupled, re- 
spectively. The new discoveries stem largely from tar- 
geted high-contrast imaging of greater numbers of YSOs 
and wide-field imaging surveys of star-forming clouds. In 
the optical and near-infrared, the broader application of 
established observing techniques, rather than any instru- 
mental improvements, has driven the recent expansion in 
the number of disks discovered. At longer mid-infrared 
wavelengths, larger telescopes equipped with more potent 
mid-infrared imagers became available. These results are 
described in this section and images of some of the new 
disks are shown in Figs. 1 and 2. A complete list of pub- 
lished YSO disks resolved in scattered light appears in Ta- 
ble 1. Most of these are around T Tau stars, although 
some are around Herbig AeBe stars. An on-line catalog 
of circumstellar disks resolved in scattered Ught, thermal 
emission, or molecular lines is now available; the URL is 
,http://www.circumstellardisks.org 

2.1. Optical and Near-Infrared Imaging 

2.1.1 Coronagraphic and Direct Imaging. The most general 
case, and also the most observationally challenging, is the 
detection of disk scattered light in the presence of direct 
starlight. Subtraction of a fiducial reference star is almost 
always necessary to reveal the disk, and a coronagraph is 



usually employed to suppress stellar diffraction. Various 
HST instruments have produced excellent images of disks 
around bright stars, while groundbased AO imaging of these 
disks is often hampered by the instability of the point spread 
function. 

The brightness of Herbig AeBe stars makes them excel- 
lent disk imaging targets, and three have been found to show 
interesting circumstellar nebulosity. HD 100546 has a disk 
viewed from high latitudes which contains what appear to 
be two spiral arms (Pantin, Waelkens, and Lagage, 2000; 
Grady et al., 2001; Augereau et al, 2001). HD 163296 
shows an inclined disk with a hint of a cleared central zone 
and bifurcation into upper and lower reflection nebulosities 
(Grady et al., 2000). The nebulosity around AB Aur shows 
a wealth of structure (Grady et al., 1999; Fukagawa et al, 

2004) , but it is unclear how much is associated with the disk 
and how much with a more extended envelope. Fukagawa 
et al (2003) report a possible disk around HD 150193A. 

The transitional disk HD 141569A has been the subject 
of a great deal of observational and theoretical work. The 
first images in the near-infrared by Weinberger et al (1999) 
and Augereau et al (1999) appeared to show a large radial 
clearing in the disk at a radius of 250 AU. Optical imag- 
ing by Mouillet et al. (2001) and Clampin et al. (2003) re- 
vealed an asymmetric spiral-like feature and showed that 
the cleared region was not completely empty. HD 14 1569 A 
is in a multiple system. Dynamical studies suggest that stel- 
lar fly-bys or a recent periastron passage by the companions 
could be the origin of the observed spiral feature (Augereau 
and Papaloizou, 2004; Ardila et al, 2005; Quillen et al, 

2005) . 

Perhaps the most significant new disk imaged around a 
directly visible star since Protostars and Planets IV is that 
of TW Hya (Krist et al, 2000; Trilling et al, 2001; Wein- 
berger et al, 2002). This face-on disk is also the closest T 
Tauri star disk, at a distance of 56 pc, and has been well- 
studied at millimeter wavelengths. The disk radial profile 
shows a sharp break in slope at a radius of 140 AU which is 
seen in multiple independent datasets. The physical origin 
of this sudden fading in the outer disk is unclear. Long-slit 
spectra of the TW Hya disk with HST STIS (Roberge et al, 
2005) find that the disk has a roughly neutral color from 
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Fig. 2. — Scattered light images of foiir edge-on or silhouette YSO disks newly resolved since 1998. Left-to-right: Edge-on disks of 
HV Tau C (Stapelfeldt et al, 2003), 2MASSI1628137-243139 {Gmsso et al, 2003), IRAS 04158+2805 {Menard et al, in preparation), 
and Orion 216-0939 {Smith et al, 2005). 



5000 A to 8000 A (see their Fig. 5), consistent with previ- 
ous broad-band imaging in the same spectral region {Krist 
et al., 2000) and near-infrared {Weinberger et al, 2002). 

Despite the interesting results highlighted above, most 
YSO disks (those whose presence is inferred from infrared 
and submm excess emission) remain undetectable in scat- 
tered light. There are probably multiple reasons for this 
fact. Disks with outer radii smaller than 0.5" , the typi- 
cal inner working angle cutoff for current coronagraphs and 
PSF subtraction techniques, would not extend far enough 
from their host star to be readily detected. This cannot be 
the whole explanation, however; the disk of MWC 480 has 
a large outer radius in milUmeterwave CO emission {Simon 
et al, 2000), but repeated attempts to detect it in scattered 
light have failed {Augereau et al., 2001; Grady et al, 2005). 
There are other examples. Many disks may simply be in- 
trinsically fainter than the ones that have been imaged in 
scattered light to date; they may be geometrically flatter 
(and thus tending toward self-shadowing), depleted of small 
particles, or have undergone chemical processing that has 
reduced their albedo. With the typical YSO disk currently 
undetectable against the direct Ught of its parent star, im- 
proved coronagraphic instrumentation will be needed to un- 
derstand the diversity in disk scattering properties and their 
underlying causes. 

2.7.2 Edge-On Disks. Edge-on, optically thick disks nat- 
urally occult their central stars and in the process present 
their vertical structure to direct view. Observations of edge- 
on disks require high spatial resolution but not high con- 
trast, so AO systems are competitive. Fiu-thermore, the ab- 
sence of stellar PSF artifacts make edge-on systems partic- 
ularly amenable to scattered light modeling. Two of the first 
known examples, HH 30 and HK Tau B, have been exten- 
sively studied over the past few years, and are discussed at 
greater length in sections 2.4, 5.1, and 5.2 below. Addi- 
tional examples are very valuable for comparative studies 
of disk scale heights, flaring, and dust properties. Finding 
new edge-on disks and imaging them across a wide range 
of wavelengths are high priorities for future research. 

The edge-on disk CRBR 2422.8 was discovered in the p 
Oph cloud core by Brandner et al. (2000a). A model for the 
source combining scattered light images and mid-infrared 



spectra is presented by Pontoppidan et al. (2005). Near 
p Oph lies another new edge-on disk, 2MASSI1628137- 
243139. Grossoetal. (2003) discovered and present models 
for this object and note a peculiar near-IR color difference 
between the two lobes of its reflection nebula. Three new 
edge-on disks have been found via adaptive optics as com- 
panions to brighter stars. Jayawardhana et al. (2002) and 
Chauvin et al. (2002) discovered LkHa 263C in a quadru- 
ple system, and Jayawardhana et al. present initial models 
suggesting a disk mass of 0.002 Mq. PDS 144 is a com- 
pact edge-on disk around a companion to an Ae star {Perrin 
et al, in preparation). The mysterious nature of HV Tau C 
was finally resolved by Monin and Bouvier (2000) to be an 
edge-on disk around the tertiary star. Modeling of HST im- 
ages of HV Tau C by Stapelfeldt et al. (2003) indicates that 
this disk also possesses a small circtmstellar envelope. Fi- 
nally, a new kind of edge-on disk has been identified around 
the Perseus source ASR 41: an extended shadow from a 
(presumably) compact disk is cast across foreground cloud 
material {Hodapp et al, 2004). 

2.1.3 Silhouette Disks. The Orion Nebula continues to 
be a unique and fertile ground for finding new resolved 
disks. Several hundred compact ionized globules ("pro- 
plyds") hkely contain circumstellar disks, but have a mor- 
phology dominated by photoevaporative processing in their 
H II region environment. Some of these contain internal 
silhouettes clearly reminiscent of disks; and there is also a 
distinct category of pure silhouette disks lacking any exter- 
nal ionization. In both cases, the disk silhouettes are visible 
in Ha images as foreground absorption to the H 11 region. 
New silhouette sources are reported in a comprehensive pa- 
per by Bally et al (2000), and by Smith et al (2005). A 
particularly interesting new source is Orion 216-0939; like 
Orion 114-426 {McCaughrean et al, 1998; Shuping et al, 
2003), this is a giant edge-on silhoutte more than 1000 AU 
in diameter, and with bipolar reflection nebulae. 

Despite the discovery of proplyds in M8 {Stecklum et 
al, 1998) and NGC 3603 {Brandner et al, 2000b) and de- 
spite HST imaging of NGC 2024, NGC 2264, M16, M17, 
and Carinae, no silhouette disks have been identified in H 
II regions other than Orion. Menard et al (in preparation) 
found the first example of a silhouette disk in nearby Tau- 



Table 1 

YSO Disks Imaged in Scattered Light 



Obiect Name 


Type 


Outer Radius 


Recent Reference 


TWHya 


face-on 


220 AU 


Roberge et al. (2005) 


HD 141569 


transition 


370 AU 


Clampin et al. (2003) 


HD 100546 


spiral 


360 AU 


Grady et al. (2005) 


HD 163296 




450 AU 


Grady et al. (2005) 


HD 150193A 




190 AU 


Fukagawa et al. (2003) 


AB Aur 


disk+envelope 


> 300 AU 


Fukagawa et al. (2004) 


CB 26 


edge-on 


380 AU 


Stecklum et al. (2004) 


CRBR 2422.8-3423 


edge-on 


105 AU 


Pontoppidan et al. (2005) 


DG Tau B 


edge-on 


270 AU 


Padgett etal. (1999) 


GGTau 


CBring 


260 AU 


Krist et al. (2005) 


GM Aur 




500 AU 


Schneider et al. (2003) 


HH 30 


edge-on 


225 AU 


Watson & Stapelfeldt (2004) 


HKTauB 


edge-on 


105 AU 


McCabe et al. (2003) 


HVTauC 


edge-on 


85 AU 


Stapelfeldt et al. (2003) 


Haro 6-5B 


edge-on 


280 AU 


Padgett etal. (1999) 


IRAS 04302+2247 


edge-on 


420 AU 


Wolf etal. (2003) 


UY Aurigae 


CBring 


2100 AU 


Potter etal. (2000) 


2MASSI J1628137-243139 


edge-on 


300 AU 


Grosso et al. (2003) 


IRAS 04325+2402 


edge-on 


30 AU 


Hartmann et al. (1999) 


OphE-MM3 


edge-on 


105 AU 


Brandner et al. (2000a) 


ASR41 


shadow 


< 3100 AU 


Hodapp et al. (2004) 


LkHa 263 C 


edge-on 


150 AU 


Jayawardhana et al. (2002) 


Orion 114-426 


edge-on 


620 AU 


Shaping et al. (2003) 


Orion 216-0939 


edge-on 


600 AU 


Smith et al. (2005) 



rus clouds: IRAS 04158+2805 shows a cone of scattered 
light, a jet, and a silhouette 3000 AU in diameter projected 
in front of diffuse Ha emission near V892 Tau. 

2.2. Mid-Infrared Imaging 

The usual objective when imaging YSOs in the mid- 
infrared is to resolve extended thermal emission from their 
inner disks. At the distances of the nearest star-forming 
clouds, extended 10-20 /xm thermal emission has only been 
detected in a handful of luminous Ae stars (HD 100546, AB 
Aur, and V892 Tau; Liu et al, 2003, 2005; Chen and Jura, 
2003; Pantin et al, 2005). In lower-luminosity T Tauri 
stars, the 10-20 /xm emssion is spatially unresolved - even 
in the relatively nearby case of TW Hya {Weinberger et al., 
2002). The absence of extended 20 fim emission in T Tauri 
stars may pose a problem for models of flared disk spectral 
energy distributions (Chiang and Goldreich, 1997). These 
models postulate a superheated disk upper layer where stel- 
lar radiation is predominantly absorbed by small dust parti- 
cles, which then radiate inefficiently at longer wavelengths. 
Large grains or flatter disks may be needed to explain the 
fact that most of these sources are unresolved in the mid- 
infrared. 

One of the more significant and surprising disk imaging 
results of recent years has been the detection of scattered 
light from three YSO disks in the mid-infrared. McCabe et 
al. (2003) found that the edge-on disk of HK Tau B appears 



as an extended source in sensitive 10 /im Keck images. The 
good alignment of this 10 fxm nebulosity with the optical 
scattered light (Stapelfeldt et al, 1998), its extent well be- 
yond a reasonable diameter for disk thermal emission, and 
its monotonically decUning flux density from 2 to 10 /xm 
argue that this is scattered light. It is not clear if the orig- 
inal source of emission is the star or the inner part of the 
disk. The fact that some edge-on disks might be seen en- 
tirely via scattered light in the mid-infrared was already in- 
dicated by Infrared Space Observatory photometry of HH 
30 (Wood et al, 2002); HK Tau B is the first resolved ex- 
ample. A second is in the case of GG Tau, where Duchene 
et ai (2004) clearly detect the circumbinary ring in deep 
3.8 jim images. By the same arguments as above, this must 
also be scattered light. A third source, PDS 144 (Perrin et 
al, in preparation), is an A star that appears as a spectacular 
bipolar nebula in 10 jim images; the relative contributions 
of scattered light and PAH emission in this image are still 
being assessed. As discussed in section 5.3 below, scattered 
Ught at these long wavelengths is a powerful diagnostic of 
large grains in circumsteUar disks. Additional examples of 
resolved mid-lR scattered light from YSO disks can be ex- 
pected in the future. 

2.3. Polarimetric Imaging 

Imaging polarimetry can confirm the presence of scat- 
tered light nebulosity and offer clues to the location of 
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embedded illuminating sources. The strength of the ob- 
served polarization depends on the dust grain properties, the 
scattering geometry, the degree of multiple scattering, and 
the polarization induced by any foreground cloud material. 
Spatially resolved polarimetry has been reported for only a 
few disks. The best example is GG Tau, where the geometry 
of the circumbinary ring is very well understood, and thus 
the polarimetric results can be readily interpreted. Silber et 
al. (2000) found that backscattered light from the far side of 
the ring was very highly polarized at 1 /xm, requiring that 
the scattered light originate from submicron-sized grains. 

A new apphcation for polarimetry has emerged in 
high contrast AO. Groundbased AO systems provide a 
diffraction-limited image core, but also an extended, uncon- 
trolled seeing halo. The hght from this halo can overwhelm 
the faint nebulosity of a circumstellar disk. Differential 
polarimetry exploits the fact that scattered light from YSO 
disks is highly polarized while the seeing halo has virtu- 
ally zero polarization. By simultaneously imaging in two 
polarizations, the unpolarized halo can be removed, and 
the polarized disk light more clearly seen. Potter (2003) 
and Perrin et al. (2004) report detections of several YSO 
nebulosities with this technique; most are circumstellar en- 
velopes. A few disks have also been studied with this tech- 
nique, notably TW Hya {Apai et al, 2004) and LkCa 15 
(Potter, 2003). The latter has not been detected in several 
conventional unpolarized imaging searches. 

2.4. Variability 

T Tauri stars commonly show photometric variability on 
timescales of a few days to months. In the youngest stars, 
these variations are now thought to be the result of hot spots, 
variation in accretion rate, and occupation by warps in the 
inner disk, all of which are natural consequences of mag- 
netospheric accretion mechanism discussed in the chapter 
by Bouvier et al. Disks also show variability, and in this 
section we describe the best studied cases, HH 30 and AA 
Tau. 

2.4.7 Photometric Variability in HH 30. Similarly to their 
central stars, disks show photometric variabiUty. The scat- 
tered light from the disk should follow as the star bright- 
ens and fades. The integrated magnitude of HH 30 varies 
on timescales of a few days over a range of more than 
one magnitude in both V and / (Wood et al, 2000). The 
range is slightly greater in V than in /. Early photometry 
suggested that the variability was periodic, but subsequent 
studies have not confirmed this. 

The range of variability of the disk is hkely to be a lower 
limit on the range of variability of the star, as the multi- 
ple optical paths taken by scattered light will likely act to 
smooth the stellar variations to some degree, both because 
of the finite speed of light (173 AU per day) and because 
the disk is illuminated by light from a range of stellar az- 
imuths. Cool spots produce stellar variability with ranges 
of no more than one magnitude {Herbst et al, 1994), and 
so cannot explain the variabiUty of the disk. Some other 
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Fig. 3. — HSTAVFPC2 images showing tlie morpliological vari- 
ability of HH 30 (Watson and Stapelfeldt, in preparation). Tlie im- 
ages are from (a) 2001 February, (b) 1995 January, and (c) 1998 
March. All images were taken in a broad filter centered at 675 nm. 
The images show the left side brighter, both sides having similar 
brighnesses, and the right side brighter. 

mechanism must be at work in HH 30, perhaps hot spots, 
variations in accretion rate, or occultations. This is consis- 
tent with the strong veiling component observed by White 
and Hillenbrand (2004). 

2.4.2 Morphological Variability in HH 30. More interesting 
and unexpected are the quantitative changes in the morphol- 
ogy of the scattered light (Burrows et al, 1996; Stapelfeldt 
et al, 1999; Cotera et al, 2001; Watson and Stapelfeldt, 
in preparation). These include changes in the contrast be- 
tween the brighter and fainter nebulae over a range of more 
than one magnitude; changes in the lateral contrast between 
the two sides of the brighter nebula over a range of more 
than one magnitude (see Fig. 3); and changes in the lateral 
contrast between the two sides of the fainter nebula over a 
range of about half a magnitude. 

The timescales for the variabihty are uncertain, but are 
less than one year. Thus, the variability is not due to 
changes in the disk structure at radii of 100 AU but rather 
changes in the pattern of illumination at radii of 1 AU or 
less. Thus, this morphological variabiUty is fascinating as 
it allows us to peek into the central part of the disk and pos- 
sibly constrain the geometry of the accretion region. 

However, before we use the variability to constrain any- 
thing, we need to understand its origin. Two mechanisms 
have been suggested; Wood and Whitney (1998) have sug- 
gested that incUned hot spots could illuminate the outer 
parts of disk like lighthouses, and Stapelfeldt et al (1999) 
have suggested that warps in the inner disk could cast shad- 
ows over the outer disk. 

In the case of HH 30, one might hope to distinguish be- 
tween the mechanisms on the basis of the color of the mor- 
phological variability. HST observations show no signifi- 
cant differences between V and /. Naively, one would ex- 
pect illumination by hot spots to produce stronger variabil- 
ity in the blue whereas optically-thick shadowing should 
be neutral in color. However, spectroscopic studies show 
that the veiling component is often relatively flat in the red 
(Basri and Batalha, 1990; White and Hillenbrand, 2004), 
so the lack of color is inconclusive. 

Determining the timescale might help to distinguish be- 
tween the mechanisms, although if the putative warp in the 
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Table 2 

Publicly Available Scattering and Thermal Equilibrium Codes 



Code 


URL 


HO-CHUNK 


|http://gemelli.colorado.edur'bwhitney/codes/codes.html 


MC3D 


|http://www.mpia-hd.mpg.de/FRINGE/SOFTWARE/mc3d/mc3d.html 


Pinball 


http://www.astrosmo.unam.mx/"' a.watson/ 


RADMC and RADICAL 


http://www.mpia-hd.mpg.de/homes/dullemon/radtrans7l 



inner disk is locked to the star as it appears to be in AA 
Tau (discussed below) this will also be inconclusive. HST 
observations show that the asymmetry can change from one 
side to the other in no more than about six months but do 
not significantly constrain a lower limit on the period or, 
indeed, demonstrate that the variability is periodic rather 
than stochastic. Since the net polarization vectors of the 
two sides of the bright nebula are not parallel (see Fig. 6 of 
Whitney and Hartmann, 1992), the asymmetry should pro- 
duce a polarimetric signature. On-going polarimetric moni- 
toring of HH 30 may help to determine the timescale of the 
asymmetry. 

2.4.3 AA Tau. AA Tau seems to be a prototype for both 
mechanisms suggested to explain the morphological vari- 
ability of HH 30, apparently possessing both inclined hot 
spots and occulting inner-disk warps {Bouvier et ai, 1999; 
Menard et ai, 2003; O'Sulivan et ai, 2005), both the result 
of an inclined stellar magnetic dipole. The disk around AA 
Tau has recently been imaged in scattered light by Grady 
et al. (in preparation), although previous observations with 
similar sensitivities did not detect the disk. It seems that 
the disk was finally detected because it was observed at an 
epoch in which the inner-disk warp was at least partially 
occulting the star. 

2.4.4 Other Disks. It would be useful to be able to compare 
the variability in HH 30 to other disks, to understand what is 
unusual and what is common. In this respect, recent second- 
epoch observations of several objects by Cotera et al. (in 
preparation) are very useful. Nevertheless, since most disks 
have been observed at only one or two epochs, there is little 
that we can say about their variability. 

3. MODELING TECHNIQUES 
3.1. Radiation Transfer 

There are several techniques for simulating scattered 
light images of disks and envelopes including the single 
scattering approximation (e.g.. Dent, 1988; Burrows et 
al, 1996; D'Alessio et al, 1999) and direct integration 
of the equation of radiation transfer under the assump- 
tion of isotropic scattering (e.g., Dullemond and Dominik, 
2004). However, by far the most common techniques are 
Monte Carlo simulations or integrations (e.g., Lefevre et 
al, 1982, 1983; Bastien and Menard, 1988, 1990; Whit- 
ney and Hartmann, 1992, 1993; Lopez, Mekarnia, and 



Lefevre, 1995; Burrows et al, 1996; Whitney et al, 1997, 
2003a,b, 2004; Lucas and Roche, 1997, 1998; Stapelfeldt 
et al, 1999; Wood et al, 1998, 1999, 2001; Lucy, 1999; 
Wolfet al, 1999, 2002, 2003; Bjorkman and Wood, 2001, 
Cotera et al, 2001; Watson and Henney, 2001; Schnei- 
der et al, 2003; Stamatellos and Whitworth, 2003, 2005; 
Watson and Stapelfeldt, 2004). Faster computers and im- 
proved algorithms allow these simulations to be fast, incor- 
porate anisotropic scattering, polarization, and fully three- 
dimensional circumstellar geometries and illuminations. 

For optical and near-infrared simulations one can nor- 
mally assume that scattered starlight dominates the im- 
ages and there is no contribution from dust reprocessing. 
With this assumption, an image at a specific wavelength 
and orientation can be calculated in a few minutes on cur- 
rent computers. Simulations at longer wavelengths must 
include thermal reprocessing and calculate thermal equi- 
librium (e.g., D'Alessio et al, 1998, 1999; Whitney et al, 
2003a,b). 

Some groups have made their scattering and thermal 
equilibrium codes publicly available (see Table 2), and 
these tools are now being used by the community. 

3.2. Density Distributions 

Modelers take different approaches to the density distri- 
bution in disks. One extreme is to assume disks are verti- 
cally isothermal, the dust is well mixed with the gas, and 
the surface density and scale heights are power laws in the 
radius. Another extreme is to solve self-consistently for 
the temperature using thermal equilibrium, to solve for the 
vertical density distribution using pressure equilibrium, to 
solve for the surface density assuming an accretion mecha- 
nism with a constant mass-transfer rate, and to include dust 
settling. In between are many intermediate approaches, for 
example, solving for thermal equilibrium and vertical pres- 
sure equilibrium but imposing a surface density law. 

The approaches are complementary. The simple power- 
law disks have very little physics, but have many "knobs" 
that can be arbitrarily adjusted to represent a wide range of 
disk density distributions and thereby cover the very real 
uncertainties in our understanding of these objects. On the 
other hand, the approaches that incorporate more and more 
physics have fewer and fewer "knobs". In one sense, they 
are more realistic, but only to the extent that our understand- 
ing of disk physics is correct. Unfortunately, there are real 
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Table 3 

Parameters Constrained by Scattered-Light Images 



Inclination Parameters 

Edge-on inclination, mass-opacity, forward scattering, and scale height 

Intermediate inclination, mass-opacity, forward scattering, and outer radius 

Face-on inclination, radial dependence of scale height 



gaps in our knowledge. For example, thermal equilibrium 
depends on the dust opacity, which is not well known; the 
details of accretion are still actively being researched, with 
disk viscosities uncertain by orders of magnitude; and disks 
may not have a constant inward mass-transfer rate. 

Both approaches are useful, but is important to under- 
stand the strengths and weaknesses of each. Simple param- 
eterized models allow one to investigate the dependence of 
other properties on the density distribution (e.g., Chauvin 
et al, 2002; Watson and Stapelfeldt, 2004), but are limited 
in what they can tell us about disk physics. On the other 
hand, models with more physics are necessary to test and 
advance our understanding of disks (e.g., D'Alessio et al., 
1998, 1999; Schneider et al, 2003; Calvet et al, 2005), 
but our knowledge of the input physics in these models is 
still incomplete. For these reasons, future modeling efforts 
will continue to tailor their approaches to modeling the den- 
sity distribution according to the specific problem being ad- 
dressed. 

3.3. Fitting 

Determining how well scattered light images can be 
modeled is important for testing physical models of flared 
disks and collapsing envelopes. Many studies have been 
succesful at reproducing the overall morphology and inten- 
sity pattem of scattered light images (e.g., Lucas and Roche, 
1997; Whitney etai, 1997; Wood et al, 2001). However, 
quantitative model fitting is more convincing and provides 
a better test of physical models. Burrows et al (1996) 
and Krist et al (2002) applied least-squares fitting tech- 
niques to single-scattering models of the disks of HH 30 
and HK Tau B and determined their physical structure and 
scattering properties. More recently, Watson and Stapelfeldt 
(2004) applied least-squares fitting techniques to multiple- 
scattering models of the disk of HH 30 and allowed the den- 
sity structure and dust scattering properties to be free pa- 
rameters. Their results provide strong constraints on the cir- 
cumstellar density, structure, and dust properties (see sec- 
tion 5.2). Glauser et al. (in preparation) have recently fitted 
models to scattered-light images, polarization images, and 
the SED of IRAS 04158-H2805. 

4. DISK STRUCTURE FROM SCATTERED LIGHT 
IMAGES 

Interpreting scattered light images of optically-thick 
disks is an inverse problem. The physics that takes us from 



a three-dimensional distribution of emissivity and opac- 
ity to a two-dimensional distribution of surface brightness 
removes a great deal of information. In this section we dis- 
cuss what information is lost and what can be recovered. 
The orientation from which a disk is viewed is a critical 
factor in this. For this reason, we classify disks as edge-on, 
intermediate-inclination, and face-on. The parameters that 
can be most directly constrained in each case are summa- 
rized in Table 3. 

4.1. Edge-On Disks 

Opaque material close to the equatorial plane of edge- 
on optically-thick disks occults the star; all that is seen 
at optical and near-infrared wavelengths are two nebulae 
formed by light scattered by material away from the equato- 
rial plane. These nebulae tend to be dominated by material 
in the segment of the disk closest to the observer. The infor- 
mation present in the images can be sunmiarized as follows: 

1. The brightness ratio of the nebulae. This is largely 
sensitive to the inclination of the disk with respect to 
the observer. When a symmetrical star-disk system is 
observed in its equatorial plane, the nebulae should 
have equal brightness. As disk is tilted with respect 
to the observer, the nebula that is closer to the Une of 
sight to the star dominates (see the top row of Fig. 4). 

2. The minimum separation of the nebulae. The min- 
imum separation occurs on the projected axis of an 
axisymmetric disk. It is largely sensitive to the to- 
tal mass-opacity product of the disk. As the mass- 
opacity product of a disk is increased, the dark lane 
separating the two nebulae widens (see the second 
row of Fig. 4). A good example of this is the de- 
creasing separation of the nebulae of HH 30 with in- 
creasing wavelength (e.g.. Figs. 6 and 8 of Cotera 
et al, 2001, and Fig. 2 of Watson and Stapelfeldt, 
2004). One cannot use the width of the dark lane to 
constrain the mass or the absolute opacity individu- 
ally without additional information. This is because 
the appearance of the disk at optical and near-infrared 
wavelengths is a scattering problem, and in the equa- 
tions governing such problems the mass density and 
opacity per unit mass always appear as a product. On 
the other hand, one can use changes in the width of 
the lane to quantify relative changes in the opacity as 
a function of wavelength (see section 5.2). 
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Fig. 4. — Scattered-light models of an edge-on disk. The center 
column shows model Al of Burrows et al. (1996). The left and 
right coltmms show models with, from top to bottom, different in- 
clinations from face-on, different mass-opacity products kM (in 
cm^ g^^ Mq, scale height normalizations Hi) (in AU), and phase 
function asyrrmietry parameters g. The contours have the same 
level in each panel and are spaced by factors of 2. See Burrows et 
al. for precise definitions of the density distribution and parame- 
ters. 



3. The increase in separation of the nebulae with in- 
creasing projected distance from the star. By this we 
refer to the degree of curvature of the boundary of the 
dark lane separating the two nebulae or the degree of 
apparent flaring of the nebulae. This is largely sensi- 
tive to the effective scale height of dust in the outer 
part of the disk (say, at the radii which dominate the 
observed scattered light). Larger scale heights lead to 
greater curvature or flaring (see the third row of Fig. 
4). Good examples of this are the relatively flared 
nebulae of HH 30, which has a scale height of 16-18 
AU at 100 AU, and the very flat nebulae of HK Tau 
B, which has a scale height of about 4 AU at 50 AU 
(see section 5.1). 

4. The vertical extent of each individual nebula. By 
this we mean the apparent height of each nebula in 
comparison to its apparent diameter. Again, this is 
largely sensitive to the effective scale height of dust 
in the outer part of the disk. Larger scale heights lead 
to more extended nebulae (see the third row of Fig. 



4). Good examples of this are again the vertically ex- 
tended nebulae of HH 30 and the narrow nebulae of 
HKTauB. 

5. The degree to which the surface brightness of the 
nebulae drops with increasing projected distance 
from the star. This is largely sensitive to the de- 
gree to which forward scattering dominates the phase 
function. Enhanced forward scattering produces a 
more centrally-concentated distribution and a more 
pronounced decrease in surface brightness with in- 
creasing projected distance (see the fourth row of 
Fig. 4). The light that emerges close to the projected 
disk axis has scattered through relatively small an- 
gles and the light that emerges far from the projected 
disk axis has scattered through relatively large an- 
gles. Increased forward scattering enhances the for- 
mer and diminishes the latter. A good example of 
this is the better fits for HH 30 obtained by Watson 
and Stapelfeldt (2004) with high values of the phase 
function asymmetry parameter g. 

However, much information that is lost or ambiguous: 

1 . The luminosity of the star. All of the light we see is 
processed through the disk, either by scattering or by 
absorption and thermal emission. 

2. The iimer region of the disk. The observed disks ap- 
pear to be sufficiently flared that the light scattered 
from the outer parts comes from a relatively small 
range of angles and passes either completely above 
the inner disk or through regions of the inner disk 
with similar extinctions. In these case, one could re- 
move the inner 30 AU of the disk and the observed 
pattern of scattered light would not change, although 
the total brightness of scattered light might. 

An exception to this general rule may be HH 30, in 
which is it is possible that parts of the inner disk shad- 
ows parts of the outer disk (see section 2.4). Another 
exception may be disks in which the inner part of the 
disk shadows the outer part completely. 

3. The outer radius. Scattered light images of edge-on 
disks sometimes show a relatively sharp outer cut- 
off and this is often interpreted as the outer radius of 
the disk. However, it is possible that the disk con- 
tinues beyond the bright nebulae, but that it is shad- 
owed. Such an extension would likely not be seen in 
scattered Ught. There is evidence that this is the case 
in the sihouette disk Orion 1 14-426, where the scat- 
tered light nebulae do not extend to the edge of the 
sihouette (see Fig. 2d of McCaughrean et al, 1998). 
An extended disk like this might also be detectable in 
CO. 

4. The radial dependence of the surface density and ef- 
fective scale height. Neither of these have strong ef- 
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1. The brightness ratio of the nebulae. This depends 
principally on the mass-opacity product and the in- 
clination (see the top two rows of Fig. 5). 




AU 

Fig. 5. — Scattered-light models of an intermediate-inclination 
disk. The center column shows model Al of Burrows et al. (1996) 
inclined to 60 degrees from face-on. The left and right columns 
show models with, from top to bottom, different inclinations from 
face-on, different mass-opacity products kM (in cm^ M©, 
and phase function asymmetry parameters g. The contours have 
the same level in each panel and are spaced by factors of 2. See 
Burrows et al. for precise definitions of the density distribution 
and parameters. 

fects on the models, at least as long as they are con- 
strained to lie within plausible ranges. Worse, even 
their subtle effects are degenerate; putting more mass 
at larger radii produces similar effects in the images 
to increasing the scale height at larger radii {Burrows 
et al, 1996; Watson and Stapelfeldt, 2004). 

It is clear from this that scattered Hght images of edge-on 
optically-thick disks provide Umited but important informa- 
tion on disk properties. Those which can be most cleanly 
separated are the inclination, the relative opacities at dif- 
ferent wavelengths, the scale height in the outer part of the 
disk, and the degree of forward scattering in the phase func- 
tion. 

4.2. Intermediate-Inclination Disks 

Intermediate-inclination disks are those in which the star 
is directly visible but both nebulae are still present (al- 
though one may be lost in the noise). In contrast to edge-on 
disks, in principle one nebulae extends over all azimuths 
from the star, although some azimuths may be lost in the 
glare of direct light from the star. Scattered-hght images 
of intermediate-angle optically-thick disks have not been as 
extensively studied as edge-on disks, but preliminary stud- 
ies {Schneider etal, 2003; Quijano, 2005) suggest that the 
following information is present: 



2. The separation of the nebulae. Like the brightness 
ratio, this depends principally on the mass-opacity 
product and the inclination (see the top two rows of 
Fig. 5). 

3. The degree to which the surface brightness drops with 
increasing projected distance from the star. This de- 
pends principally on the the degree to which forward 
scattering dominates the phase function (see the third 
row of Fig. 5). 

4. The ratio of scattered hght to unscattered light. That 
is, the relative brightness of the nebulae and the star. 
This depends on the scale height in the outer disk, 
the albedo, the degree to which forward scattering 
dominates the phase function, the inclination, and the 
mass-opacity product. (A complication here is that 
light scattered very close to the star is difficult to dis- 
tinguish from direct light from the star.) 

5. The outer radius. The presense of two nebulae in 
many systems suggests that there is not a cold, col- 
lapsed, opticaUy-thick disk that extends significantly 
beyond the bright nebulae. This allows one to con- 
strain the outer radius of the optically-thin disk with 
some degree of confidence. For example, GM Aur 
shows two nebulae which suggest that the true outer 
radius of the optically-thick disk is around 300 AU 
{Schneider et al., 2003). 

Again, much information is lost: 

1. The inner region of the disk. In current images, in- 
formation on the inner disk is lost under the PSF of 
the star, even when observing in polarized light (e.g., 
Apai et al, 2004). Infrared interferometry, covered 
in the chapter by Millan-Gabet et al, can provide a 
great deal of information on the very innermost part 
of the disk at radii of less than 1 AU. Coronagraphic 
images from space or, perhaps, from extreme AO 
systems will be required to recover information on 
scales larger than those available to interferometers 
but smaller than those lost under current PSFs. 

2. The radial dependence of the surface density and ef- 
fective scale height. As with edge-on disks, neither 
of these have strong effects on the observed morphol- 
ogy. 

Scattered-light images of intermediate-angle optically- 
thick disks are thus more difficult to interpret than those of 
edge-on disks. The only parameter that seems to be cleanly 
separated is the outer radius of the disk. The mass-opacity 
product and the inclination are to a large degree degenerate 
(compare, for example, the first two rows of Fig. 5, which 
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Fig. 6. — Scattered-light models of a face-on disk. The center 
column shows model Al of Burrows et al. (1996) inclined to 10 
degrees from face-on. The left and right columns show models 
with, from top to bottom, different inclinations from face-on and 
different indicies /3 in the scale-height power law. The contours 
have the same level in each panel and are spaced by factors of 2. 
See Burrows et al. for precise definitions of the density distribution 
and parameters. 



show similar changes in the brighter nebula), but if millime- 
ter data or an SED are available to constrain the inclination, 
the mass opacity product can be obtained from scattered- 
light images. The dust asymmetry parameter can be con- 
strained if disk structure is known with adequate certainty. 

4.3. Face-On Disks 

Face-on disks are those in which departures from radial 
symmetry are dominated by inclination (combined with the 
phase function) rather than vertical structure (at least for ax- 
isymmetric star-disk systems). The only example we have 
of such a disk around a YSO is TW Hya {Krist et al, 2000; 
Trilling et al., 2001; Weinberger et al., 2002; Apai et al., 
2004; Roberge et al, 2005). 

The information present in scattered-light images of 
face-on optically-thick disks is: 

1. The ellipticity. This is by definition determined by 
the inclination and phase function in axisymmet- 
ric star-disk systems. Departures from constant el- 
lipticity can provide a fascinating insight into non- 
axisymmetric illumination or disk structure. For ex- 
ample, Roberge et al (2005) show evidence that the 
disk around TW Hya is significantly elliptical be- 
tween 65 and 140 AU and essentially circular beyond 
this region; they suggest that this suggests that the 
disk is warped. (See the first row of Fig. 6.) 

2. The dependence of the surface brightness on radius. 
This depends essentially on the radial dependence of 
the scale height. If the disk is well-mixed, this pro- 
vides information on the radial dependence of the 
temperature. (See the second row of Fig. 6.) 



3. The ratio of scattered light to unscattered light. That 
is, the relative brightness of the nebula and the star. 
This depends principally on the scale height in the 
outer disk, the albedo, the phase ftmction, and the 
mass-opacity product. 

Information on the inner disk is not present for the same 
reasons as in the case of intermediate-inclination disks. In- 
formation on the outer disk radius is not present for the 
same reasons as in the case of edge-on disks. 

5. HIGHLIGHTS OF APPLIED DISK MODELING 

In this section, we present four examples of scientific 
results from the modeling of scattered-light disks around 
YSOs. Our first three are related to the evolution of dust, 
which is also discussed in the chapters by Natta et al and 
Dominik et al. The fourth concems evidence for an inner 
hole in the disk around GM Aur, which may have been 
cleared by a planet (see the chapter by Papaloizou et al). 
While many disks have been imaged, the cases discussed 
below are perhaps the best examples of how disk structure 
and dust properties have been derived from modeling of 
multi-wavelength image datasets. 

5.1. Evidence for Dust Settling 

The gas scale height in disks is determined by the bal- 
ance between the stellar gravitational force which com- 
presses the disk toward its midplane and gas pressure 
which acts to puff the disk up. In vertical hydrostatic 
equilibrium in a vertically isothermal disk, the gas scale 
height is related to the local gas temperature according to 
H{r) = .JkT{r)r^/GM^m, where m is the mean molec- 
ular weight of the disk gas. Direct determination of the 
disk temperature is difficult, but reasonable results can be 
obtained by modeling the infrared and millimeter spectrum. 
A complication is that disks likely have vertical temperature 
inversions (Calvet et al, 1991, 1992; Chiang and Goldre- 
ich, 1997; D'Alessio etal, 1998). 

If the dust and gas are well mixed, then the dust follows 
the gas density distribution and will have the same scale 
height. As we have seen in section 4.1, the dust scale height 
can be determined from high-resolution scattered-light im- 
ages. It is important to note that the height of the scatter- 
ing surface above the disk plane is not equal to the dust 
scale height; rather, the scattering surface defines the locus 
of points with optical depth unity between the star and in- 
dividual disk volume elements, and can be any number of 
dust scale heights above the midplane. The primary obser- 
vational indicators of the scale height are the vertical extent 
of the nebula, the degree of curvature of the nebulae, and the 
sharpness of the dark lane. Disks with large scale heights 
appear more vertically extended ("fluffy") and have curved 
disks whereas those with small scale heights are vertically 
narrow and have almost parallel nebulae. Exact value for 
the scale height must be determined by fitting models to the 
images. 
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Table 4 
Scale Heights 



Object 


H at 50 AU 


^equivalent 


Reference 


HH 30 


6.3 AU 


51 K 


Watson and Stapelfeldt (2004) 


IRAS 04302+2247 


6.1 AU 


48 K 


Wolfetal. (2003) 


HV Tau C 


6.5 AU 


35 K 


Stapelfeldt et al. (2003) 


GM Aur 


3.4 AU 


18 K 


Schneider et al. (2003) 


HKTauB 


3.8 AU 


8K 


Stapelfeldt et al. (1998) 



The comparison of the scale heights for the gas and dust 
offer a unique opportunity to test the assumption that disks 
are vertically well mixed. Burrows et al. (1996) derived an 
equivalent temperature from the dust scale height for the 
HH 30 disk at 100 AU radius under the assumption that the 
disk was vertically well mixed. The result was broadly con- 
sistent with expectations from a simple thermal model, and 
suggested that the assumption that the disk was vertically 
well mixed was correct. 

Since Protostar and Planets IV dust scale heights have 
been derived for several other YSO disks. In two cases, HV 
Tau C, and IRAS 04302+2247, initial scale height deriva- 
tions implied unreasonably high equivalent temperatures 
and it is difficult to imagine why the dust would be more 
extended than the gas {Stapelfeldt et al., 2003; Wolf et al, 
2003). These appear to be caused by the presence of cir- 
cumsteUar envelopes in addition to the disks, with the en- 
velopes producing more diffuse nebulosity than expected 
for a pure disk. By adding an envelope to the density 
distribution, it was possible to remove this effect, and de- 
rive dust scale heights more consistent with simple thermal 
models. Scale height values for these and three other disks 
are shown in Table 4. The scale heights in the original ref- 
erences have been extrapolated to a reference radius of 50 
AU to facilitate comparison. 

The dust scale heights measured to date fall into two 
groups with the values differing by almost a factor of two. 
The three objects with larger scale heights all have out- 
flows indicating ongoing accretion, whereas the two ob- 
jects with smaUer scale heights have Uttle or no outflow 
activity. This could be observational evidence that accret- 
ing disks are systematically more puffed-up (warmer) than 
non-accreting disks. The smaU equivalent temperature im- 
plied for the HK Tau B disk merits particular attention. Un- 
fortunately the infrared spectral energy distribution for this 
disk is incomplete, so this value cannot be compared to a 
well-constrained disk thermal model. However, it is highly 
unlikely that the gas in the disk could actually be as cold as 
8 K. Instead, this could be a case where the assumption that 
the dust and gas are well mixed is not correct. Instead, it ap- 
pears more likely that the dust has decoupled from the gas 
and partially settled toward the disk midplane, an expected 
stage of disk evolution {Dubrulle et al. 1995; Dullemond 
and Dominik, 2004). In this case, the equivalent tempera- 



ture derived from the dust scale height would be lower than 
the true gas temperature. 

The vertical scale height of dust in a YSO circumstellar 
disk could be a key indicator of its structure and evolution- 
ary state. It would be very valuable to accumulate scale 
height measurements for larger numbers of edge-on disks, 
to see how unique the HK Tau B results are and uncover the 
full diversity of scale heights in the YSO disk population. 

5.2. Evidence for Dust Growth from the Opacity Law 

In optically thin media, the color of scattered light di- 
rectly depends on the wavelength dependence of the grain 
opacity. Small grains like those in the ISM have a much 
higher scattering cross section at shorter wavelengths, and 
thus optically thin nebulae will have strongly blue colors 
relative to their illuminating star. However, YSO disks are 
optically thick in the optical and near-infrared. In this situ- 
ation the color of the scattered light no longer depends on 
the grain opacity, and the disk will appear spectrally neutral 
relative to the star - even if smaU grains are the dominant 
scatterers. Neutral colors have been observed for the disks 
of TW Hya (Krist et al, 2000; Weinberger et al, 2002; 
Roberge et al, 2005) and GM Aur {Schneider etal, 2003), 
consistent with this expectation. A mechanism that can pro- 
duce non-stellar colors in reflection from an optically thick 
disk is a wavelength-dependent dust albedo or phase func- 
tion, but this is thought to be a small effect. Given these 
considerations, how can scattered light from an optically 
thick disk be used to constrain the grain properties? 

The answer is to study changes in nebular spatial struc- 
ture as a function of wavelength. For small ISM-like parti- 
cles, the T — \ scattering surface is located in lower-density 
regions above the disk midplane at optical wavelengths, and 
shifts into higher-density regions nearer the midplane at 
near-infrared wavelengths. Conversely, large grains acting 
as grey scatterers would produce a reflection nebula whose 
spatial structure would not vary with wavelength. 

For edge-on disks, the key observable is a narrowing of 
the central dust lane as the object is imaged at progres- 
sively longer wavelengths (see the second row of Fig. 4). 
This behavior is seen in Orion 114-426 (McCaughrean et 
al, 1998), IRAS 04302+2247 {Padgett et al, 1999), HH 
30 {Cotera et al, 2001), and HV Tau C {Stapelfeldt et al, 
2003). It provides clear evidence that the scattering dust 
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Fig. 7. — A comparison of data and models for the GG Tau cir- 
cumbinary ring, in the / band (0.8 yUm; top row) and M band (3.8 
(Um; bottom row). The ring is inclinced 37° from face on, with the 
forward scattering edge projected above the binary, and backscat- 
tering side projected below the binary. The key observable is the 
brightness contrast between the front and back sides. The first 
model has a maximum grain size of 0.3 /xm. It matches the / band 
image well, but in the M band it predicts too much backscatter- 
ing and not enough forward scattering. The second model has a 
maximum grain size of 0.9 fim. In this model, the backscattered 
flux becomes too small in the / band, and still too large in the 
M band. A reasonable mach to the M band image is obtained in 
the third model (maximtmi grain size of 1.5 /xm), but the I band 
backscattering is underpredicted. No single grain size distribution 
accounts for the phase function effects at both wavelengths. Fig- 
ure and results from Duchene et al. (2004). 



grains are dominated by small particles. But how small? 
A quantitative answer can be derived by fitting scattered 
light models to multiwavelength image data sets. This was 
first done in the case of HH 30 by Cotera et al. (2001). 
From modehng of HST near-infrared images, they found 
that the dust lane thickness changed less quickly with wave- 
length than expected for standard interstellar grains, and 
interpreted this as evidence for grain growth in the disk. 
Watson and Stapelfeldt (2004) extended this analysis by in- 
cluding optical images and by considering a wider range of 
possible density structures, but found essentially the same 
result: the most likely ratio of grain opacities between 0.45 
and 2.0 fim is 2.0 for the HH 30 disk, versus a value of 10 
expected for ISM grain models. 

While the disk of HH 30 appears to show some grain 
evolution, modeling of multiwavelength images of the 
edge-on disks of HV Tau C (Stapelfeldt et al, 2003) and 
IRAS 04302H-2247 (Wolf et al, 2003) finds grain opac- 
ity ratios consistent with standard ISM grains. Both of 
these sources possess circumstellar envelopes in addition 
to disks, so the presence of primitive grains could reflect 
this ongoing infall from the ISM onto the disks. The disk 
of HK Tau B shows only subtle changes in its dust lane 
thickness between optical and near-IR images, and may 
represent a more evolved system; a firm conclusion on its 
dust properties can be expected from future model fitting. 
Dust properties in the giant edge-on silhouette disk Orion 
114-426 are uncertain; Throop et al. (2001) found the ra- 
dial extent of the silhouette to be achromatic between 0.66 
and 1.87 /xm; Shaping et al. (2003) found that it was chro- 



matic between 1.87 and 4.05 iim; while McCaughrean et 
al. (1998) showed that the dust lane thickness between the 
lobes of reflected light was clearly chromatic between 1.1 
and 2.0 /xm. Additional modeling of this source is needed. 
Observations and modeling of a broader sample of edge- 
on disks offer an opportunity to probe the diversity of dust 
properties across the variables of disk age, disk enviroment, 
and accretion signatures, and should be vigorously pursued. 

5.3. Evidence for Dust Stratification from tlie Pliase 
Function 

GG Tau is a binary T Tauri star with 0.3" (42 AU) pro- 
jected separation. It hosts the most prominent example of a 
circumbinary ring of dust and gas. The ring has been stud- 
ied in the millimeter lines and continuum (Guilloteau et al., 
1999; Woodet al., 1999), near-infrared scattered hght (Rod- 
dier et al., 1996; Wood et al, 1999; Silber et al, 2000; Mc- 
Cabe et al, 2002), and optical scattered light {Krist et al, 
2002, 2005). Through these studies, the density structure of 
the GG Tau ring is perhaps now the best understood of all 
YSO disks. A key feature of the ring is its intermediate in- 
clination of 37° from face on. This spatially separates both 
the foreground and background parts of the ring from each 
other and from the central binary. This "clean" configura- 
tion allows the relative strength of forward scattering and 
backscattering to be directly measured. This quantity can 
be a powerful diagnostic of dust properties in the ring. 

The scattering phase function strongly favors forward 
scattering when the grain size is comparable to the wave- 
length. At wavelengths much larger than the grain size, 
scattering becomes more isotropic. A comprehensive study 
of phase function effects in the GG Tau ring was recently 
carried out by Duchene et al (2004). Using new images 
taken at 3.8 pra, the longest wavelength to date at which 
scattered hght from the ring has been detected, and existing 
images at shorter wavelengths, DucMne et al modeled the 
wavelength dependence of the phase function from 0.8 to 
3.8/um. Highlights from their results are shown in Fig. 7. 
The results show that the 3.8 ^,m scattered light must arise 
from dust grains larger than those in the ISM, whereas the 
0.8 /im scattered light must simultaneously originate from 
much smaller particles. Additional evidence for small par- 
ticles is provided by the 1.0 fxm polarimetry of Silber et al 
(2000), who found the backscattered light was highly polar- 
ized and thus dominated by submicron grains. 

The GG Tau phase fiinction results indicate that no sin- 
gle power-law distribution of grain sizes can simultaneously 
account for the observations. Duchene et al (2004) suggest 
that these results can be explained by a vertically stratified 
disk in which the grain size increases towards the midplane. 
In this view, large grains responsible for the 3.8 /um scat- 
tered light would be located in a denser region closer to 
the disk midplane, a region that shorter wavelength photons 
cannot reach. The stratification could be due to dust settling 
to the ring midplane. Alternatively, it might be the case that 
the preference for scattering at grain sizes comparable to 
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the wavelength is so strong that the large grains dominate 
the 3.8 fim scattered hght, even though they are less numer- 
ous than the small grains. In that case, vertical dust setthng 
would not necessarily have taken place. 

This picture is very attractive. However, Krist et al. 
(2005) point out that current models of GG Tau do not fuUy 
explain the observations. For example, no current model 
simultaneously reproduces the total brightness of the disk, 
its color, and its azimuthal variation. Furthermore, the ra- 
tio of brightness of the near and far parts of the disk has 
been observed to vary with time. Additional observations 
and models will be required to confirm the suggestion of a 
stratified disk. It would be very valuable to perform a sim- 
ilar study in other disk systems; unfortunately, the GG Tau 
ring is currently unique. 

5.4. Combining Scattered Light Images and SEDs 

GM Aur has a disk viewed at an intermediate inclination. 
The disk is clearly seen in WFPC2 and NICMOS images af- 
ter careful subtraction of a reference PSF. Scattered light is 
detected between 0.4" and 2.1" (55 and 300 AU) from the 
star. Modeling of the scattered light images by Schneider et 
al. (2003) derived a dust vertical scale height of 8 AU at a 
radius of 100 AU and a disk inclination 56° from face-on, 
and demonstrated that there was also a renmant circumstel- 
lar envelope. However, the coronagraphic occulting spot 
used for the observations blocked any information on the 
properties of the inner disk. To access that region, Schnei- 
der et al. turned to the infrared spectral energy distribution 
(SED). Emission between 2 and 70 /xm probes the disk tem- 
perature structure inside 40 AU, and thus is complementary 
to the results of the scattered light modeling. Using the disk 
model derived from the scattered light images, and appro- 
priate assumptions about grain properties, Schneider et al. 
were able to reproduce GM Aur's SED, including the mil- 
limeter continuum points (see Fig. 8). The very small near- 
IR excess emission requires that the irmer part of the disk 
be optically thin. Rice et al. (2003) showed that the inner 
region may have been cleared by a Jovian-mass planet. 

More recently, Calvet et al. (2005) have presented 
Spitzer IRS spectroscopy of GM Aur. This data requires 
an inner optically-thin region that extends to 24 AU from 
the star. This is much larger than previous estimates based 
on modeling of broad-band mid-infrared photometry. The 
Spitzer IRS data show that the inner disk is not empty, but 
contains a small amount of small dust grains which produce 
the silicate emission feature. To fit the long wavelength 
SED, the average grain size in the outer disk must be larger 
than in the inner, optically-thin region. 

Scattered light images do not generally provide informa- 
tion on the inner part of the disk. In addition to IRAS data, 
high-quality mid-infrared spectra from the Spitzer Space 
Telescope are becoming available for many YSO disks. 
Scattered light images reduce many of the degrees of free- 
dom in modehng disk spectra and SEDs, leading directly to 
more robust SED models. Future combined studies of mul- 




Wavelength (|im) 

Fig. 8. — Scattered light images and SEDs of the classical T 
Tauri star GM Aur (Schneider et al., 2003). The top left panel 
shows the scattered light image and the top left a model. The 
lower panel shows the SED. The dark line is a model for the stellar 
photosphere emission. The white line shows a model combining 
the stellar emission and the disk excess emission, with disk model 
parameters taken from scattered light results. For the photometry 
references, see Schneider et al. (2003). 

tiwavelength disk images, continuum spectra, and SEDs 
will provide critical tests for models of disk evolution. 

6. PROSPECTS FOR FUTURE ADVANCES 

6.1. Science Goals 

The science goals for scattered-light imaging of disks 
around YSOs center on gaining a better understanding of 
disk structure, disk evolution, and dust evolution in the con- 
text of planet formation. 

A key goal is to determine an evolutionary sequence for 
dust, especially dust growth and changes in shape and com- 
position. One manifestation of dust evolution is to change 
the dust opacity law. Work in HV Tau C and HH 30 has 
shown that the opacity law can be determined from HST im- 
ages in the optical and near-infrared using current modeling 
codes. Extending this work is simply a matter of obtaining 
the relevant observations and "cranking the handle". Obser- 
vations of thermal emission at longer wavelengths will help 
to restrict the degeneracies in the scattered light modeling. 
We would expect to find disks with different opacity laws, 
some with near-ISM opacity laws like HV Tau C and others 
with less opacity laws hke HH 30. The science will come 
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from comparing the opacity laws with models for dust evo- 
lution and from correlating the opacity laws with system 
properties such as age, accretion activity, disk mass, and 
stellar multiplicity. 

The evidence for dust settling discussed in section 5.1 
is as yet tentative. To advance, we need to be able to con- 
fidently compare the equivalent temperatures derived from 
the dust scale heights to the real temperatures in the disk or, 
more generally, the vertical distributions of the dust and the 
gas. This will require detailed calculations of the structure 
of the disk that solve for the vertical structure under the as- 
sumption of thermal equiUbrium but allow for uncertainties 
in the radial structure. Again, including mid-infrared and 
millimeter data will be very useful. 

Of course, we are ultimately interested in detecting 
young planets in disks, and especially those planets that 
are still accreting. Direct detection will be difficult (see the 
chapters by Beuzit et al. and Beichman et al), but we can 
hope to observe the gap cleared by a planet as it grows. At 
the moment, these are most clearly seen in infrared spectra, 
and the contribution of scattered fight images is mainly to 
reduce the number of degrees of freedom in the disk geom- 
etry and inclination and thereby improve the refiabiUty of 
the models for the infrared emission. 

We would like to understand the transition between 
optically-thick protoplanetary disks and optically-thin de- 
bris disks in order to determine the timescale over which 
planetesimals and gas - the building blocks of planet for- 
mation - are present. Advances are being made in this 
field with Spitzer imaging and spectroscopy. Scattered light 
imaging can contribute because small amounts of dust can 
potentially be detected even at radii at which infrared detec- 
tions are difficult, although the realization of this goal will 
require advances in high-contrast imaging. 

6.2. Instrument Advances 

The biggest instrumental advance for the study of YSO 
disks will be the advent of the Atacama Large Millimeter 
Array (ALMA). With spatial resolution surpassing that of 
HST, very high resolution spectroscopy for chemical and 
kinematic studies, and the sensitivity to study the disks of 
low-mass stars as far away as Orion, ALMA will have a ma- 
jor impact on the field. However, even in the era of ALMA, 
scattered light imaging will still make significant contribu- 
tions to disk studies. First, scattered fight traces the surface 
where stellar photons deUver energy to the disk. It will still 
be necessary to characterize this interaction region if the 
disk temperature structure, and thus its chemical nature, is 
to be understood. Secondly, while millimeter wavelengths 
probe large particle sizes in the dust population, scattered 
light imaging provides information on the small particles. 
Both are needed to provide a full picture of disk grain prop- 
erties and their time evolution. The combination of mil- 
limeter maps and scattered light images at comparable reso- 
lutions will be a powerful synthesis for disk science efforts. 

The imaging performance of large groundbased tele- 



scopes can be expected to continue its evolution. Partic- 
ularly important will be new extreme adaptive optics in- 
struments at the Gemini and VLT observatories. Their im- 
proved contrast performance should enable additional de- 
tections of YSO disks in the near- and mid-infrared. The 
maturation and increasing application of differential imag- 
ing polarimetry should also yield exciting results. Several 
concepts for extremely large (D ~ 30 m) telescopes are now 
being studied. When realized (maybe not until after Proto- 
stars and Planets VI), these faciUties will provide a three- 
fold advance in spatial resolution. Higher resolution images 
will improve our knowledge of all aspects of disk structure. 
The inner holes in systems such as TW Hya and perhaps 
GM Aur should be resolvable. An exciting possibility is 
the detection of hot young planets near to and perhaps even 
within YSO disks and characterization of their dynamical 
interactions. 

Among future space missions in design and develop- 
ment, two will provide important capabilities for disk scat- 
tered light imaging. The NASA/ESA James Webb Space 
Telescope will have a 6.5 m primary mirror and operate 
from 1 to 28 /xm. It will be a superb telescope for imag- 
ing disks in mid-infrared scattered light and will provide 
roughly 0.2" resolution. In the near-infrared, it will offer 
almost three times better resolution than HST, but is un- 
likely to provide improved contrast. From the point of view 
of high contrast imaging, the Terrestrial Planet Finder Coro- 
nagraph mission would be extremely exciting, as it would 
be able to detect scattered light from disks as tenuous as 
our own solar system's zodiacal light. However, as this 
misssion may be more than a decade away, several groups 
are proposing smaller coronagraphic space telescopes that 
might be realized sooner. 

While waiting for these future developments, the ex- 
panded application of existing scattered light imaging capa- 
bifities (AO and mid-IR imaging from large ground-based 
telescopes and high-contrast imaging with HST) should 
continue unabated. 

6.3. Modeling Advances 

The new observations described in this review and the 
observation advances outfines above suggest a wealth of 
detailed data on circumstellar disks will become available. 
What will this data demand from the codes and models? 

Future codes will need to produce high-resolution im- 
ages and integrated spectra at wavelengths stretching from 
the optical to the millimeter They should be able to model 
three-dimensional distributions of sources and opacity and 
should incorporate accurate dust scattering phase functions, 
polarization, and the effects of aligned grains. Evidence 
for dust growth and sedimentation requires that codes no 
longer restrict themselves to homogenoeus dust properties, 
but must be able to treat multiple dust species with different 
spatial distributions. They will probably solve for radia- 
tion transfer using Monte Carlo techniques in optically thin 
regions and the diffusion approximation in optically thick 
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regions. Computers are expected to become increasingly 
parallel in the future. Many codes can already run in par- 
allel, but those that cannot will need to be modified to do 
so. In this context, Monte Carlo algorithms have the ad- 
vantage over classical algorithms as they often have natural 
parallelism. Many radiation transfer codes now have the ca- 
pability to determine the density structure from given disk 
physics or incorporate density structures from dynamical 
simulations. These will be increasingly useful in providing 
detailed tests of disk structure. On the other hand, parame- 
terized density models will continue to be a useful tool for 
mitigating our incomplete knowledge of disk physics. 

The combination of future data and fiiture codes will al- 
low us to study dust properties, dust settling, disk struc- 
ture, disk-planet interactions, accretion, and disk evolution. 
In the longer term, advances in techniques for numerical 
simulations coupled with increases in computing power and 
parallel processing make accurate radiation hydrodynamic 
simulations of disk formation and evolution a distant but 
realistic goal. 

Perhaps the most important near-term work is to apply 
scattered light modeling techniques across already extant 
disk image datasets. Observers have provided a significant 
number of new, high-quality images, but the corresponding 
modeling efforts have not kept pace. 
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ported by the Centro de Radioastronomia y Astrofisica of 
the Universidad Nacional Autonoma de Mexico, HST GO 
program 9424 funding to the Jet Propulsion Laboratory of 
the California Insitute of Technology, and the Programme 
National de Physique Stellaire (PNPS) of CNRS/INSU, 
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